Finite size effects in Barabasi- Albert growing networks 



O 

o 

(N 

X) 
D 

\o 

(N 



I 

o 
o 



(N 
> 
00 
(N 

ON 

o 

\o 
o 



X3 
O 

o 



X 



B. WaclawH and I. M. Sokolov 
Marian Smoluchowski Institute of Physics, Jagellonian University, Reymonta 4, 30-059 Krakow, Poland, 
Institut fur Physik, Humboldt Universitdt zu Berlin, Newtonstrasse 15, 12489 Berlin, Germany 

We investigate the influence of tfie network's size on tfie degree distribution tt^ in Barabasi- Albert 
model of growing network with initial attractiveness. Our approach based on moments of -Kk allows 
to treat analytically several variants of the model and to calculate the cut-off function giving finite 
size corrections to Hk. We study the effect of initial configuration as well as of addition of more 
than one link per time step. The results indicate that asymptotic properties of the cut-off depend 
only on the exponent 7 in the power-law describing the tail of the degree distribution. The method 
presented is this paper is very general and can be applied to other growing networks. 



I. INTRODUCTION 



Complex networks have been widely studied by physi- 
cists and researchers in many other fields. It was pointed 
out that many natural and social networks are scale-free 
(SF), i.e. that the degree distribution Hk giving the prob- 
ability that randomly chosen node has exactly k nearest 
neighbors behaves as ~ for sufficiently large degrees 
fc, typically with the exponent 7 laying between 2 and 
4 [H 0]- This most striking feature of real- world com- 
plex networks is well described within theoretical models 
proposed during the last decade. However, in any finite 
network the power-law behavior of the degree distribu- 
tion TTfc can hold only for the values of k small enough 
in comparison to the number of nodes N . Both experi- 
mental data and theoretical models of scale- free networks 
indicate that the behavior of ivk for fc ^ 1 for a finite net- 
work exhibits two regimes: below some kmax it follows 
the power-law behavior as in an infinite network while 
above k^ax it displays a much faster decay. The charac- 
teristic degree k^-^ax which separates these two regimes is 
called a cut-off. Intuitively, the cut-off is due to the fact 
that since the overall number of links present in a finite, 
non-degenerated graph is bounded from above, so is also 
the degree of each node. The truncation of the power- 
law affects many properties of networks, especially the 
percolation ones, important in models for infection prop- 
agation of real diseases or computer viruses. We note 
here that the cut-offs due to finite size effects have to 
be distinguished from the distributions' changes when a 
small sarnple of the network is considered, as discussed 
in Ref. Q, the fact which makes the analysis of real 
networks very involved. 

Many attempts were undertaken to estimate the po- 
sition of the cut-off for different network models. Here 
we consider sparse networks with fixed average degree 
(which corresponds to 7 > 2), and focus on the cut-off 
stemming from finite-size effects. They depend strongly 
on whether the network is growing and thus exhibits cor- 
relations age-degree, or is a homogeneous one Q, ob- 
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tained in a sort of thermalization process leading to sta- 
tistical equivalence of all nodes. The results for homoge- 
neous networks [1, 0, 0] indicate that kmax should scale 
as ~ Ari/(7-i) for 7 > 3 and as - TV^/^ for 2 < 7 < 3 
as a result of structural constraints leading also to the 
occurrence of correlations between node degrees. These 
results based either on probabilistic arguments or on ex- 
treme values statistics did not provide an explicit form 
of the degree distribution for a finite network. On the 
other hand, the statistical ensemble approach of Ref. Q 
lead to kmax scaling rather like N'^/^^^^) for 2 < 7 < 3, 
at variance with the prior results. 

The models of growing networks allow for a more rig- 
orous treatment. In Q a simple model of growing tree 
network has been solved exactly. The authors have de- 
duced the form of corrections to the degree distribution 
and found that the cut-off scales as ~ N'^/h-'^) for any 
exponent 7 larger than two. The scaling for 7 < 3 is 
thus different from that observed for homogeneous net- 
works and can be understood as a result of degree-degree 
correlations [6|. 

In the present paper we give a general approach to ob- 
taining degree distributions in growing networks of finite 
size. It is based on moments of connectivity distribution 
TT/j and works for any case for which one is able to write 
the rate equation for tt^ depending on the size N . Be- 
cause of its generality, the method is interesting by itself. 
In the paper we discuss its application to the Barabasi- 
Albert (BA) model with initial attractiveness ao, intro- 
duced in different contexts several times in the past [ld| . 
The model is defined as follows. Starting from a finite 
graph of tlq nodes, at each step a new node is introduced 
and joined to m previously existing nodes with probabil- 
ity being proportional to fc -I- ag, where k is the degree of 
a node to which new link is established. This rule, called 
preferential attachment, causes that nodes with high de- 
grees get connected to more and more nodes. The sim- 
plest version of that model with ao = was discussed 
by Barabasi and Albert 11]. An effective algorithm of 
building such a network with the attachment probability 
proportional to k consists in focusing on the links, and 
not on the nodes of the graph. Each link ij is then con- 
sidered as a couple of oriented links i j and j i 
pointing into opposite directions. Because the probabil- 
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ity of picking at random an oriented link pointing onto 
the node with degree k is proportional to k, the preferen- 
tial attachment rule is simply realized by choosing that 
node as the one to which the newly introduced node is 
attached. 

The model has been extended to the case of ag 7^ and 
solved in the thermodynamical limit in [T^ . The popu- 
larity of this model is partially explained by its three 
important properties: i) the generated network is scale- 
free with tunable exponent 7 = aq -I- 3, ii) the generated 
network is a small world one, iii) the algorithm can be 
easily implemented on a computer. The preferential at- 
tachment rule can be generalized to non-linear attach- 
ment kernels [isl . [l^ . These however do not give pure 
scale-free behavior and thus will not be discussed here. 

Some of the results we present here were obtained pre- 
viously in Ref. l9|. Our approach is however different 
from the one of [9| and allows for solving more sophis- 
ticated variants of the model. The method is explained 
in Section II where it is introduced on the example of 
the simplest BA model with m — 1 and ao — 0. In 
Section III we consider the influence of the initial con- 
figuration on the cut-off function. This is of particular 
importance since it has been pointed out 1151 1 that the 
functional form of the cut-off depends strongly on the 
seed graph. We show however that the leading behavior 
of the cut-off function is independent on initial conditions 
and for large k the cut-off can be well approximated by 
a Gaussian form. Section IV is devoted to generalized 
attachment kernel k + ag. Now both the exponent 7 7^ 3 
and the cut-off depend on Uq. We calculate the cut-off 
function explicitly and show that it is no longer Gaussian 
but exhibits the scaling predicted in 0. Section V deals 
with non-tree graphs. Some mathematical formulas we 
exploit in the paper are discussed in Appendices. 



II. FINITE SIZE EFFECTS IN A GROWING 
SCALE-FREE TREE 



For the presentation of our method we begin with the 
simplest variant of the model, namely the BA tree net- 
work with m = I, no — 2 and without initial attrac- 
tiveness. The mathematical approach developed in this 
Section will be applied to more involved situations in 
Sections III and IV. 

We start from putting down the rate equation for 
Nk{N), the average number of nodes with degree k for 
the network consisting of N nodes. For simplicity we 
assume that the growth process starts from a dimer con- 
figuration, that is TiQ = 2 nodes joined by a link. Thus, 
iVfc(2) = 2Sk,i, where 5kq stands for the Kronecker delta. 
At each time step a new node is added and connected 
to a node i chosen at random out of N old vertices with 
probability proportional to its degree fci, thus obtaining 
a new network with -I- 1 nodes. The network is a tree 



graph. The recursion formula for Nk{N) has the form 



+ 



Nk{N + 1) ^ Nk{N) + 6ks + 
fc — 1 k 



NkiN). (1) 



The equation is exact for any N. The Kronecker delta 
stands for the introduction of a new node with degree 
one. The next term corresponds to the appearance of a 
node with degree fc as a result of link's addition to the 
node with degree k—1. The last term describes how fast 
nodes with degree k disappear due to link's addition to 
those nodes. 

In the limit A — > 00 the solution is given by A^^ w 
Ntt^^ where, according to Ref. p^ . 



k{k + l){k + 2) 



(2) 



In this paper we are interested in the general solution for 
Nk{N) = Nnk{N), with TTk{N) being the degree distri- 
bution of a finite network. Our method to obtain this dis- 
tribution is based on using its moments. We shall do this 
in several steps. First, we consider moments for 7rfc(A^), 
which are of particular importance for many problems 
including epidemic spreading and percolation 16!|. Let 
Mn{N) be the n-th moment of A^fc(A^) and Hn{N) be the 
n-th moment of the degree distribution Trk{N): 



Mn{N) = ^fc"A^fe(iV) 



fe=i 



Y,k^MN)^-rMn{N). 
fe=l 



1 

'n' 



(3) 
(4) 



Multiplying both sides of Eq. ([T]) by /c" and summing 
over k we get now recursive relations between moments 
Mn{N) for different A^. For instance, the equations for 
the first three moments read: 



A/o(^+ 1) 
Mi{N+l) 

M2{N+l) 



Afo(A) + 1, 
2N -I 

2(A^^^^(^) + 1' 

2A 
2(A- 1) 



M2{N) 



Ml (TV) 
2(A- 1) 



(5) 
(6) 

1,(7) 



with the initial condition Af„(2) = 2 imposed by Eq. ^ 
for the dimer configuration. The important feature of 
these (and higher) relations is that for any n, Mn{N + 1) 
depends only on lower moments. We can thus solve these 
equations for Mn{N) by forward substitution, starting 
from the lowest moment. Equation (O has an obvious 
solution Mq{N) = N. Solution of Eq. ([6]) is also easy 
to guess, being just twice the number of links: Mi{N) = 
2{N — 1). In order to solve Eq. ([7|) we insert the solution 
for Ml and we solve the resulting equation by standard 
methods. We get 



M2iN) = 2{N - 1)H{N - 1), 



(8) 
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where H{n) = X]r=i 1/* is the harmonic number. Higher 
moments can be computed using similar tricks, but the 
formulas get much more complicated. The appropriate 
moments for TTk{N) are 



Mo = 1, 

III = 2-2/N, 

fi2 = {2-2/N)H{N 



(9) 
(10) 
(11) 



We see, that ^2 ~ 2 log diverges for large N, but is still 
quite small for networks of order 10^ nodes, the largest 
real-world structures like the WWW. This suggests that 
real SF growing networks exhibit strong finite size effects 
which is indeed true jTi|| and the assumption that 112 = 00 
for the BA model is far from reality even for very large 
networks. 

The equation for fi2 gives a rough estimate for the po- 
sition of the cut-off in the degree distribution as imposed 
by finite size effects. Let us assume the cut-off in the 
degree distribution to be sharp enough and take tt/j to be 
TTk — T^k^ for k smaller than kmax and vanishing above 
kmax- Then the second moment behaves as Alogkmax- 
Comparing this with exact Eq. ([TT|) we get kmax ~ V^- 
The position of the cut-off scales therefore as a square 
root of the network size which in fact has been pointed 
out by several authors 0, 0| . We will show further explic- 
itly that such scaling indeed holds for BA model without 
initial attractiveness. 

To find the correction to ttJ^^ in the finite network 
we introduce a function Vk{N) defined as Vk{N) = 
Nk{N)/iT^^ . The function Vk{N) gives finite size cor- 
rections. From the asymptotic properties of iik in the 
limit of large we know that Vk {N) cx N for k ^ kmax 
and falls to zero far above the critical kmax- Equation 
(fT|) can be rewritten in terms of Vk{N): 



Vk{N) 



2 ■ 2{N -2) ' 
4:-2N + k , 



(12) 



To proceed, we define moments m„(-/V) for the distribu- 
tion Vk{N): 



m„(iV) = 



1 



N 

k=l 



(13) 



with the value of the normalization constant l/(iV — 1) 
motivated by later convenience. Multiplying now both 
sides of Eq. (|12p by fc" and summing over fc = 1 , . . . , 00 
we get 

1 / 

m„(Ar+l) = ^ 3 + ^c„m,(7V)-h 



+ {2N + n + l)mn{N) ] , 



(14) 



where the initial condition 

m„(2) = 3 



(15) 



stems from the configuration of the starting graph and 
the coefficients c„i are given by 

c„o = 3, and c„i = ^(^^ (i " l) ^'^^ * ^ ^' ^^^^ 

By solving equations for the first few moments we can 
infer the general solution: 



n+l 



Br, 



r(7V) ^ r(2 -f z/2) 



T{N + i/2), (17) 



where the coefficients Bni still have to be found and the 
factor r(2 + i/2) is singled out for convenience. Inserting 
this to Eq. (|14p we get after some manipulations: 



n+l 



i=0 



T{2 + i/2) 

n 71—1 



^C„,S,or(iV)+^ Cn,B,, 



i=l j=i-l 



TiN + t/2) 

r(2 + V2) ■ 



(18) 



Comparing terms of the same order in F-functions we get 
two recursion relations for Bni: 



1 / ""^ 



3=0 



for i — and 



B„ 



^ n — l 



(19) 



(20) 



for < i < n. The third recursion relation which com- 
pletes the set comes from the initial condition P^ : 



-^n,n+l — 3 ^ ^ Bjii^ 



(21) 



i=0 



and Bqq = —3. Returning to mn{N) wc find that for 
large N the leading term reads: 

mn{N) ^ i?n,n+i ^^^^^^ ^'nii^ ^ N^A,,, (22) 

with An = B„.n+i/T (■^-^)- Such scaling indicates that 
for large N, 



Wfc(iV) ~iVw(fc/V7V), (23) 
where w{x) is a function having An as its moments: 



/>oo 

An = / dxw{x)x^ . 
Jo 



(24) 
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Therefore the degree distribution for large but finite BA 
tree network is well approximated by 



(25) 



The functional form w{x) of the cut-off can be found 
analytically by reconstructing it from its moments An, 
which express through the coefficients B„.n+i- However, 
in order to calculate them one has to know all Bni result- 
ing from Eqs. (dH) , (gD]) and Fortunately, they can 
be conveniently stored by means of appropriately defined 
generating functions, for which we can write differential 
equations and solve them. In Appendix A we show that 



An — 



{2 + nynl 
m + n)/2y 



(26) 



The behavior of the cut-off function ■w{x) for large val- 
ues of the argument is dominated by the behavior of An 
for large n which can be easily obtained using Stirling's 
formula for the factorial and for the Euler's Gamma func- 
tion: 



logAn « -n\ogn. 



(27) 



Let us now compare it with moments /„ of the function 
exp[-(x/a)"]: 

x" exp [-(a;/fT)"] dx = T —— 

a \ a y 

For large n the leading term of log /„ w (n log n) /a is the 
same as in Eq. (|27|) with a — 2, which means that the 
tail ofw{x) falls like a Gaussian function. The parameter 
a can be established by comparing sub-leading terms in 
/„ and An ■ The value a = 2 obtained in this way will be 
confirmed below by direct calculation of w{x). 

From the values of the coefficients An, the function 
w{x) can be restored using a trick involving inverse 
Laplace transform. In Appendix B we show how to get 
w{x) expressed as an infinite series: 



w(x) 



n=l 



,2n+l 



(-l)"n2 



i!22"(2r 



1 



(29) 



One can check that this result corresponds to a Taylor se- 
ries for an expression involving the complementary error 
function erfc(z): 



w[x) — erfc(x/2) H ^_ e 



V4 



47r 



(30) 



the resul t gi ven in [9j, which is close to an approximate 
result of [T^ . The series in Eq. is rapidly convergent 
and, if truncated at some rimaxi can serve for numerical 
calculations. 



III. ARBITRARY INITIAL CONDITION 

In this section we generalize the method to the case 
when the network's initial configuration is a complete 
graph with more than two nodes. We proceed exactly as 
before. First we write the rate equation for the average 
number of nodes of a given degree. Second, we factor- 
ize the solution into two terms: the one corresponding 
to the thermodynamical limit N oo and another one 
Vk{N) giving the correction for finite N. Then we find 
moments for Vk{N) and deduce the form of w{x) from 
their asymptotic properties using formulas given in Ap- 
pendices. Although we still add one link per one added 
node, we shall consider here a more general rate equation 
for Nk{N): 



NkiN + 1) 



Nk{N) 
k 

2{N -I 



5k'. 



k-l 
2{N -u) 



(31) 



with two free parameters m and u. The parameter m 
corresponds to the number of links added with one new 
node. In the cases considered in the previous and in the 
present section we have m = 1. The case m > 1 is con- 
sidered in Sec. [V] The parameter u in the denominators 
in Eq. (j3ip is introduced to guarantee the normalization 
of the probability of preferential attachment and depends 
on the form of the preferential attachment kernel and on 
initial conditions. For the linear attachment kernel k the 
normalization condition is ^ kNk{N) = 2L, with L be- 
ing the overall number of bonds. We define uj to fulfill 
the equation: 2L = 2{N — lo). If we start with a complete 
graph with no nodes then lo = no(3— 7to)/2, which in par- 
ticular gives 2L = 2{N — 1) for the dimer-configuration 
from Sec. II. Here we first consider Eq. (j3ip for arbitrary 
m and uj and put appropriate values of these parameters 
at the end of computations. 

In the thermodynamical limit the degree distribution 
is given by [Tsj : 



7roo(fc) 



2m{m + 1) 
k{k + l){k + 2)' 



(32) 



for all fc > m and does not depend on lo. Like in the 
previous section, we insert Nk{N) — nj^^Vk{N) and get 
the equation for Vk{N). Then, defining the moments 
mn{N) = k'^Vk{N)/ {N — lu) we obtain the recursion 
relation: 



m„(A^+ 1) 



m„(no) 



n-1 



m]m 



c„imi{N)' 



i=0 



+ (2N - 2uj + 3 + n)mn{N) ] , 

(»o - l)"o 
2m(m + l)(no — lo) 



(33) 

(no - 1)", (34) 



where the initial condition (|34p has now a more compli- 
cated form than before. The general solution for mn{N) 
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can still be written by means of Gamma functions: 



mn{N) = 



1 



T{N + l-uj + i/2) 



T(no + l-^ + V2)' 
(35) 

which imposes that Vk{N) scales as a function of a; = 
k/^/N exactly like in Eq. (|23p. The equations for Bni 
look similar to Eqs. (Ull), and but with the 
constant 3 changed for (2 + m)m'^T{nQ + 1 — w) in the 
first and for m„(no)r(no + 1 — w) in the third equation. 
Proceeding as in Section II we find generating functions 
for Bni and finally, the expression for An'- 



A-n 



r(l + TlQ — u) 

r(no + 3/2-w + n/2) 
TOo(Tio)r(2 + no + n 
r(no + 2) 



r(3 + m 



(n + l)r(m + 2) 



(36) 



The generating function M(z) will therefore be a sum 
of two /-functions defined in Eq. (jBSp . see Appendix 
B. Using the results of Appendix C we can immediately 
write the expression for the cut-off function w{x): 



w{x) — T{l + no—Lu) 
mo{no) 



r(m + 2) ■^^^"r) + i-'-'.3-l-m,2(^) 



r(nn + 2) 



,l(^) 



(37) 



where /-functions are defined in Appendix C and lo, mo 
have been given above and depend only on uq. In figure[T] 
we plot w{x) calculated from Eq. ([57| together with the 
results of numerical simulation for finite size networks. 
One readily infers a strong dependence of w{x) on the 
size of the seed graph no. This sensitivity to the initial 
conditions has just been reported in as well as in Ref. 
(T7j where another quantity was measured. However, if 
one compares the asymptotic behavior of A„'s: 



log A„ 



1 



-nlogn 



(l-log2). 



(38) 



with that of Eq. ([281) . one immediately finds that for 
large n the function w{x) still behaves like exp(— a;^/4), 
independently of no. Therefore the degree distribu- 
tion for the BA tree model without initial attractiveness 
has always a Gaussian cut-off whose position scales as 

~ iVl/2. 



IV. GENERALIZED ATTACHMENT KERNEL 

The considerations above can be easily extended to the 
case of preferential attachment kernel k + ao with initial 
attractiveness ao > —1. In Ref. it is shown that the 
degree distribution in the thermodynamical limit reads: 



^ (2 + ao)r(3 + 2ao) r{k + ao) 



r(l + ao) r(/c + 3-H2ao) 



(39) 




FIG. 1: The cut-off function wlx) giving finite-size corrections 
to the degree distribution, calculated from Eq. (|37p for the 
starting graph with no = 3 (black lower line) and 5 nodes 
(black upper line) . Our analytical results agree very well with 
'w{x) obtained from averaged degree distributions for 2 x 10* 
generated networks of size iV = 10* (gray lines). 



For large fc, tt^ decays according to a power-law ~ k~^ 
with 7 = 3-|-ao. Therefore this model can reproduce any 
exponent 7 observed in real-world networks. In [T3| it has 
been shown that the model is equivalent to the growing 
network with re-direction (GNR), which is a convenient 
method for simulating the model on a computer. At each 
step we choose at random one node i from the set of ex- 
isting nodes. A newly introduced node is attached with 
probability 1 — r to this i, and with probability r to its 
ancestor, i.e. to the node the proper, oriented link of i 
was attached to when this node was itself introduced into 
the network. With the choice r = l/(ao -I- 2) the GNR 
model is equivalent to the BA tree model with initial at- 
tractiveness. In all numerical simulations showed in this 
section the GNR model is used. A generalized version of 
the GNR model which works for m > 1 links added per 
one new node is described in flS*] . 

The recursion formula for Nk{N) is slightly modified 
in comparison to Eq. ([1]) and reads: 



Nk{N + l) 



NkiN) 



k + ao — I 
N{2 + ao)-2 



Nk-i{N) 



ao 



-NkiN)+Sk,i, 



N{2 + ao)-2"^'" ' ' "^•'^ 

but the initial conditions remain the same as in Section 
II. Inserting now Nk{N) = Vk{N)n'^ we get the equation 
for the function Vk{N) and then the equation for the 
moments m„(A^): 

Tn„{N + 1) = h 7 r— X 

Ar(2+ao) [{N-l){2+ao)+ao) N 

Cmm,{N) + {N{2 + ao) + ao + 1 + n)m„(7V) 



i=0 



(41) 
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where c„i's read now: 

Cno = 3 + 2ao, 
Cm = (3 + 2ao) 



1 



, for i > 1. (42) 



Following the same steps as in Section II we can find 
the general functional form of mn{N) and its leading 
behavior for large networks: 



m„(iV) 









r 




"P ( 1 ao+n+l 1 
^ 2+ao ) 


{N-l) 



O _i_ QQ+n+l 
^ 2+ao 



(43) 



Therefore the func- 



with An = i3„,„+i/r 
tion Vk{N) obeys the following scaling rule: 

Vk{N) Nw {k/N^' \ , (44) 



where the function w{x) has moments An depending on 
flQ. Equation (|44l) indicates that the cut-off scales as 
where 7 = 3 + ag is the exponent in the power- 
law for vr^ . This is in agreement with arguments given in 
0. For a given size N , the cut-off shifts to lower values of 
TV when the exponent 7 increases. This implies that the 
power-law in the degree distribution can hardly be seen 
for 7 > 4, because even for large networks with = 10^ 
nodes the cut-off corresponds to the value of k^ax ^ 100 
and the power-law lasts only for 1 — 2 decades in k. This 
partially explains the fact that the SF networks with 7 
above 4 are practically unknown [2] . 

As before, we can find a general formula for An and 
estimate its leading behavior: 



log An 



1 + ap 

2 -f ao 



n log n. 



(45) 



Comparing this to Eq. ([28)1 as it has been done before, 
one obtains that for large x the function w{x) decays like 
exp [— (a;/cr)"] with 



7 ■ 



1 + ao 7 - 



1 



(46) 



This agrees very well with numerical findings and means 
that the cut-off for 7 3 is no longer of a Gaussian type. 
For 2 < 7 < 4, as often found in real networks, a is always 
larger than 1.5 and therefore the cut-off w{x) cannot be 
approximated by a pure exponential decay observed in 
some networks [2J; exponential cut-offs most likely have 
a different origin than the finite size effects. The value 
of a can also be obtained from subleading terms and is 
given by 



(2 + ao)(l + ao)-i/". 




FIG. 2: Plots of w{x) calculated from Eqs. (gH]) and (fB23)) for 
ao = —0.5, and 1 (solid lines from top to bottom) which cor- 
respond to 7 = 2.5, 3 and 4, respectively. The curves become 
flat with increasing ao. The thick gray lines are w{x) obtained 
from averaged degree distributions for 2 x 10* networks of size 
iV = 10''. The tails decay like exp(-2;°) with a = 3,2,3/2, 
respectively, in agreement with numerical findings. 



The result for M{z) has now the form of Eq. (jB5|) with 



a 



2r(i + 2ft) 

r(5 + 2ao) 
1 



ao 



ao 



6 + 3ao, b= 12 + 13ao-l-4ag, 
2, (48) 



,/3=?^,e = 4 + 2ao,C 



(47) 



and therefore w{x) is given by Eq. (|B23p with the cor- 
responding values of parameters. In Fig. [2] we plot 'w{x) 
for ao = —1/2, and 1. For the purpose of numerical cal- 
culations all series have been truncated, with the error 
being less than 10~^ in the plotted area. These results 
show that the curves become more flat with increasing 
ao and agree well with w{x) obtained in simulations of 
finite size networks. 

As before, the starting graph has large influence on the 
exact form of w(a;). We do not consider here the depen- 
dence on the size no of the seed graph, however, just like 
in Sec. Ill, one can show that asymptotic properties of 
the cut-off function are insensitive to no and therefore for 
X being sufficiently large, w{x) ~ exp [— (x/cr)"] depends 
only on ao i.e. only on the exponent 7 in the power-law 

TTfc ~ k~^. 



V. BA MODEL WITH m > 1 - BEYOND THE 
TREE MODEL 

In the previous sections we have considered growing 
networks restricting ourselves to the case when graphs are 
essentially trees (the possible cycles can only stem from 
the seed graph). The general case of rn > 1 links joined 
to a new node in a single step is much more complicated. 
Each of m proper links of a newly introduced node has to 
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be connected to one of A'' preexisting nodes according to 
the preferential attachment rule. However, since multiple 
links are not allowed, the nodes to which links have been 
connected on this step have to be excluded from the set 
of nodes available for further linking. Thus, when a new 
link is introduced, the probabilities of attaching it to one 
of the preexisting nodes depend whether the link is the 
first, second, etc. of m. The rate equation for Nk{N) 
can still be obtained in this case. However, its structure 
is very sophisticated and it is highly non-linear in N^, 
and k. This makes impossible the immediate application 
of our method. To proceed, we note that for large N 
and m <^ N the probability of choosing the same node 
as a candidate for link's attachment more than once in 
a single time step is very small. We can thus write the 
rate equation fully disregarding the exclusion of multiple 
connections, since their relative number should be small 
in the thermodynamic limit N —^ oo. For instance, for 
m = 2 the approximate equation for Ni.{N) reads: 



NkiN + 1) = Nk{N) + 5k,2 + 
NkiN). 



1 



27V- 3 



Nk-iiN)- 



2N-3 



(49) 



We assume additionally A^fc(3) = ■i6k,2 that is we start 
from a complete graph with uq — 3 (triangle). In the 
general case m > 2 the approximate equation takes the 
form of Eq. ([?T|) . The origin of all terms is the following. 
The Kronecker delta gives the addition of a node with 
m links at each time step. The number of links L = 
m{N — Lu) with 



uj — no{2m + 1 — nQ)/{2m) 



(50) 



gives the normalization factor 2L for the probability of 
attachment. The factor m hidden in 2L cancels out with 
m possibilities of choosing links at each step. For m = 2 
and no = 3 it reduces to Eq. (|49|) . 

Like before, we expect some dependence on starting 
graph, but as long as asymptotic properties of ^(a;) are 
concerned, its particular choice does not seem to be im- 
portant. It is therefore tempting to assume no = m + 1 
because it allows to have tt^ = for all fc < m at each 
time step. However, we must remember that the Eq. 
pip with Lv given by Eq. (|50p is a reasonable approxima- 
tion of the unknown true rate equation for Nk{N) only 
ii m <^ N a.t each stage of the network's growth. Thus 
our expressions for w{x) works well only for m <ti hq. 

We can check the validity of Eq. ((3T|l calculating the 
moments /i„(A^) of the degree distribution TTfe. The equa- 
tion (|3ip has a form similar to Eq. ([T]) and thus all results 
are similar to those presented in Sections II-IV. However 
the approximation we made includes only leading terms 
in the limit of large N and the quality of that approxima- 
tion can only be checked by computer simulations. For 
instance, for uq = m + 1 the first three moments /i„ (N) 



N 


fi2 exact 


fi2 from Eq. (|54)| 


error in% 


200 


34.16 


35.26 


3.2 


1000 


43.73 


45.15 


3.2 


2000 


47.83 


49.34 


3.1 



TABLE L Comparison between exact and approximate value 
of the second moment fi2 in case of m = 2. 



are 



Mo = 1, 

/ii = m{2- {m + l)/N) ^2L/N, 

m{m + 1) 



M2 



m+1 




7V-1 

E 



2i-m+l 



(51) 
(52) 

(53) 



First two moments are exact although we have used the 
approximate equation. For large N the last expression 
reduces to 



m+l 
N 



H ^-^ ( 7 - log 2- 



h(^) + log(2iV-m- 1) 



(54) 



where H{rL) are harmonic numbers and 7 ~ 0.5772 is 
the Euler-Mascheroni constant. Thus the second moment 
grows like m{m-\- 1) log A'^ and similarly to what has been 
done in Sec. II one can infer that the cut-off scales like 
\/N . In table [T] we compare the quality of formula ([M]) 
for m — 2 with fi2 found for numerical solutions of exact 
equation for Nk{N). The good agreement confirms that 
the approximation we have made performs well. 

Let us go now to the corrections to the degree distri- 
bution. In the thermodynamic limit, tt^ is given by Eq. 
(15^ . The degree distribution for large but finite networks 
corresponds to tt^ multiplied by w{k/\/N). From Sec. 
HI we know that w{x) is expressed by Eq. (|37p where m 
and rto are now completely arbitrary and w is given by 
Eq. ((50)) . The asymptotic form of w{x) ^ exp(— x^/4) 
resulting from analytical formula, is the same as in Sec. 
III. In figure [3] we compare our approximate analytical 
solution with results of direct numerical simulations for 
different no. A small deviation between analytical and 
numerical curves is evident. This deviation is the largest 
for no = 3 and the smallest for ng — 15, confirming that 
the performance of our approximation is the better, the 
larger is the seed graph. 



VI. CONCLUSIONS 

We have investigated finite size corrections to the de- 
gree distribution tt^ in the Barabasi- Albert model of 
growing network, assuming preferential attachment ker- 
nel fc -|- Oo. Imposing that tt^ is a product of known 
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FIG. 3; The function w{x) for m — 2 and for no — 3, 7, 10 and 
15 (curves from the flattest to the most peaked). The agree- 
ment between the analytical solution and simulation results is 
not as good as before due to approximate character of solved 
equation. The simulation data were obtained for A'^ = 4000. 
The dip at about x = 1 and the peak at about a; = 2 es- 
pecially pronounced in the case of larger seed graphs mean 
that for large starting graph much more nodes with high de- 
grees is present than it would be expected for the asymptotic 
power-law behavior of . 



APPENDIX A 



Let us define the following generating functions for 



n=2 



and 



(Al) 



(A2) 



for i > 0. Note that Eq. (|A1|) is not a special case of 
Eq. We start from Eq. multiplying its both 

sides by {n + l)x"/n! and summing over n — 2, . . . , oo. 
As a result one can easily show that ^0(2:) satisfies the 
differential equation: 

(e"= - l)F[,{x) = (2 - 3e"=)Fo(x) - 3{3x - l)(e"= - 1) - 3x, 

(A3) 

with Fo(0) = 0. Equation (jA3| has the following solu- 
tion: 



Fa{x) = 2-3a;-e" 



-2x _ ^-x^ 



(A4) 



distribution tt^ ^ k for infinite network and the func- 
tion w{x) giving the finite size cut-off in the rescaled 
variable x = k/N 1-'^ , we have reduced the problem to 
finding the moments of w{x). We have shown that for 
sufficiently large N they grow faster than exponentially 
and give the explicit formula for their leading behavior. 
We have found the expression for w{x) and argued that 
far above the cut-off it behaves as ^ exp(— (a;/cr)") with 
a = (7 — 1) 7(7— 2) for any exponent 7 lying in the range 
2 < 7 < 00. We discuss the dependence of w{x) on 
the initial configuration of the network and show that al- 
though it has strong infiuence on exact form of w{x), its 
asymptotic properties are independent on the seed graph. 
Moreover, we give a general expression for 'w{x) in terms 
of convergent power series in x and present plots for some 
special cases. 

The approach we have developed here allows for treat- 
ing the general BA model in a new fashion. This ap- 
proach can be also applied to other models of growing 
networks provided the recursion equations for the con- 
nectivity Nk{N) are known. 
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In a similar way, from Eq. (|20[) and using the definition 
of Cni from Eq. (|16p . we get the differential equation for 
F,{x): 



{l-e^)Fi{x) = (3e--2-z)f^,(x)-Hi?.-i,.— — X 

(i - 1)! 



3e"-3+^ i(e"-l-a;) 



(A5) 



This complicated equation has a simple solution due to 
the fact that Fi [x) contains only powers of x higher than 
i - 1: 

F^{x) = B,^,/- V4t^ — (A6) 



Now, if we multiply Eq. (j2ip by and sum over 

71 = 1, ... , 00, we find: 



Y,'^B^,n+i+Y,F,{x)^Z{e'' -x-l). (A7) 



i=0 



Substituting Fi{x) by Eqs. (|A4p and (|A6p and defining a 
new generating function for the coefficients Bn^n+i- 
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2: 

^ — ^ n ' 



(AS) 



we get the expression: 
e-^^'Gil - e-^) = 3e^ + l + e'^^ + e"^ - 66"^^, (A9) 
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as follows from Eq. (IA7[) . This equation for G{z) can be 
readily solved: 

G{z) = 3(1 -z)-^+{l- + (1 - + (1 - - 6. 

(AlO) 

The values of Bn^n+i can be obtained by using Cauchy's 
formula as integrals over a contour encircling the point 
z = in the complex plane: 

nl f dz , (2 + n)2(3 + n) 

Bn,n+^ = 2^ / ^^(^) - 2 ■ 

(All) 

Finally, from the definition of An = n+i/F (^^) we 
get Eq. 



APPENDIX B 

To get the functional form of the cut-off w{x) we define 
a generating function: 



M{z) = Y,A^- 

n=0 

Comparing this definition with Eq. 
M{z) = exp{zx)w{x)dx so that 



(Bl) 
we see that 

M{-z) = / exp{-zx)w{x) dx (B2) 
Jo 

is the Laplace transform of w(a:). Therefore w{x) is given 
by the inverse Laplace transform of M{z) or equivalently 
by the Fourier transform of M{—iz): 



w{x) 



1 

2^ 



dzp>"'A.f(-iz). 



Using the explicit form of coefficients A„ we get 

(2 + n)F(n + 3) ^„ 



M{z) = J2 



F(n + 2)F((3 + n)/2) 



(B3) 



(B4) 



This series has an infinite radius of convergence. The 
function M{z) given by Eq. l|B4p is a special case of 
more general power series: 



M{z)^J^Y. 



F(n + C)F(an + /3)' 



(B5) 



which depends on six parameters a, b, a, (3, i^, C and where 
M gives the appropriate normalization. We introduce an 
auxiliary function fa.i3,s,,Q{z)'- 



F(n + 



„^^F(n-|-C)F(an + /3) 



(B6) 



This one is a special case of a Fox- Wright Psi function 
^9*1 or of a Fox H function (see e.g. [lol). Its Fourier 
transform can be written as follows: 



dze''''f{-iz) = Ic{x)-Is{x), 

> 

(B7) 



where Ic{x) and Is{x) are cosine and sine transforms of 
real and imaginary part of f{—iz), respectively: 



1 f°° 

Is{x) = — / sin{zx)'^f{~iz)dz 
Jo 



Ic{x) 



cos{zx)^f{—iz)dz. 



(B8) 
(B9) 



We consider here only Ic{x) since repeating all calcu- 
lations presented below for Is{x) one can check that 
Ic{x) — —Is{x)- The function Ic{x) can be expressed 
via the Mellin transform f^{s) of 'Rf{~iz): 

1 1 

Ic{x) ^—7^ /5(s)r(l - s)sm{^s/2)x'ds, 
TTX 2m J^,^ 

(BIO) 



where 

fc'is) = / z'-^nf{-tz)dz. 
Jo 

The real part of f{—iz) reads: 

(-z2)«F(2n + ^) 



n=0 



r(2n + C)T(2an + (3) ' 



(Bll) 



(B12) 



We now substitute two of three Gamma functions by 
their integral representations: 



/•OO 

r(2n + = / e"*t2n+«-irfi^ 
Jo 



r(2n + c) 



(B13) 
(B14) 



where "H" denotes the Hankel contour. Changing vari- 
ables z— >w:z = 'i;j, we get the following integral 
representation: 



/■oc /• 

/5(s)=/ dte~H^''-^ 

Jo J Ha 

/•oc 

X / dww"-!^ 



due^u^ ^ X 



2^rl 



-gr(2an + /3)- 



(B15) 



The last integral is now just the Mellin transform of the 
Mittag-Leffler function defined as 



^ — ^ jTL 

^ r(2Q:n + p) 

n—0 ' 



From [2l| p. 301 we have 

„ . 2. Mellin r(s/2)r(i - s/2) 



2r(/3- 



(B16) 



(B17) 



and after simple algebraic manipulations we finally arrive 
at: 



r(e - 5)r(g/2)r(i - s/2) 
2r(c-s)r(/3-as) ■ 



(B18) 
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Now we are able to calculate f{x) = 2Ic{x) as an in- 
verse Mellin transform. Using the identity for Gamma 
functions 



r(t)r(i-t) = 



sin(7ri) ' 

and applying the residue theorem we get: 



(B19) 



r(^-^)r(i-.) 

T{(i-as)T{Q-s) 



E 



ress 



r(^-s)r(i-s)_^,_i 



Lr(c-s)r(/3-as)' 



(B20) 



where the sum runs over all points s„ at which either 
r(l — s) or r(^ — s) has a pole. The above formula sim- 
plifies for ^, C, being positive integers m, k: 



constant term vanishes under the transform, so that we 
finally get: 

w{x) = TV {axfa,f3-a,e,-i,c-i{x) + - a)/Q,/3,?,c(^)) ' 

(B23) 

where /'s are given either by Eq. (|B2l|) or by more gen- 
eral Eq. (|B20)) . 

Applying now Eq. (|B23|1 to TV = 1, a = 1, 6 = 2, a = 
1/2,^ 3/2,^ = 3,C = 2 as stems from Eq. (HH, we 
get the function w{x) for the BA tree model: 



w{x) = a;/i/2,i,24(a;) + /i/2,3/2,3,2(a;) 



E 

n=0 



{-xY 



l-n 



r(l/2-n/2) r(l-n/2) 



(B24) 



n=0 



(m — 2 — n)(7n— 3 — n)---(/e— 1— n) 
r(/3 — Q — an)n! 



(B21) 

where we make use of the fact that reSz=_„ T{z) — 
(— l)"/rt!. Let us come back to M{z) which can be ex- 
pressed through f(z) and its derivative f'{z): 

M{z) = M (az/;^^^,^(z) -f 6/.,/5,c,c(^)) • (B22) 

Then w{x) is given by Fourier transforms of the two func- 
tions /(z) and zf'{z). Integrating by parts we can change 
from ^ f Q, 

,/3,«,c(^) *° /a,/3-a,4-i,c-i(2;) +const, where the 



With help of Eq. (|B19|) we get after some algebraic ma- 
nipulations: 

/ \ 1 ^ 2n+l/ INK 2-r('^+l/2) /T3r,C^ 

w{x) = 1 > X ^ ( — 1) n — — . (B25) 

^ ' TT ^ ^ ' T(2n + 2) ^ ' 

Finally, having in mind that 



r(2n + 2) = ^n!(2n-^-l)r(n-f 1/2), (B26) 



we get Eq. 
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